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Abstract. U(l) gauge fields are decomposed into a monopole and pho- 
ton part across the phase transition from the confinement to the Coulomb 
phase. We analyze the leading Lyapunov exponents of such gauge field 
configurations on the lattice which are initialized by quantum Monte Carlo 
simulations. We observe a strong relation between the sizes of the monopole 
density and the Lyapunov exponent. Evidence is found that monopole fields 
stay chaotic in the continuum whereas the photon fields are regular. 



1. Motivation 

A main question in QCD is the existence of special classes of gauge fields 
which are responsible for confinement. In the framework of the dual super- 
conductor, color magnetic monopoles arise as the appropriate candidates. 
Their investigation necessitates a projection of the non-Abelian gauge the- 
ory on the compact U(l) theory. In lattice computations it was demon- 
strated that monopoles account for more than 90 percent of the string 
tension in the confinement and they drop toward zero across the phase 
transition. 

The study of chaotic dynamics of classical field configurations in field 
theory is motivated by phenomenological applications as well as by the 
understanding of basic principles. The role of chaotic field dynamics for 
the confinement of quarks is a longstanding question. Here, we analyze the 
leading Lyapunov exponents of compact U(l) configurations on the lattice. 
The real-time evolution of the classical field equations was initialized from 
Euclidean equilibrium configurations created by quantum Monte Carlo sim- 
ulations. The U(l) gauge field is decomposed into a monopole and photon 
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part. This way we expect to learn some details between the appearance of 
monopoles and the strength of chaotic behavior in lattice simulations. 

2. Monopoles in compact quantum field theories 




We have investigated 4(i U(l) gauge theory described by the action 

S{Ui} = J2i^-cos9p), (1) 
p 

with Ui = Ux,^, = ex.p{iex,n) and Op = Bx,,j. + 0x+fi,u - Ox+o,^ - Ox,u {v [x] ■ 
At critical coupling /3c~1.01U(l) gauge theory undergoes a phase transi- 
tion between a confinement phase with mass gap and monopole excitations 
for P < Pc and the Coulomb phase with a massless photon for (3 > Pc, see 
Fig. 1. 

We are interested in the relationship between monopoles of U(l) gauge 
theory and classical chaos across this phase transition. Following Refs. [1, 
2, 3, 4], we have factorized our gauge configurations into monopole and 
photon fields. The U(l) plaquette angles Ox^nu are decomposed into the 
"physical" electromagnetic flux through the plaquette Ox,ii.v and a number 
mx,ij,v of Dirac strings passing through the plaquette 

^x,ixv — (^x,ixv "I" 27r TRx^jxv , (2) 

where 9x,fiu £ (— +7r] and rrix^^u 7^ is called a Dirac plaquette. Monopole 
and photon fields are then defined in the following way 

x' 
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^SIJT* - - X) <^a;,a;' d'u Ox',v,x ■ (4) 
x' 

Here & acts on x', the quantities are defined in Eq. (2) 

and Ga; is the lattice Coulomb propagator. One can show that Ox,^ = 
0™°° + is up to a gauge transformation identical with the original 6x,ij, 
defined by Ux,i^ = exp(i6'a;,^). 

3. Classical chaotic dynamics from quantum Monte Carlo initial 
states 

Chaotic dynamics in general is characterized by the spectrum of Lyapunov 
exponents. These exponents, if they are positive, reflect an exponential di- 
vergence of initially adjacent configurations. In case of symmetries inherent 
in the Hamiltonian of the system there are corresponding zero values of 
these exponents. Finally negative exponents belong to irrelevant directions 
in the phase space: perturbation components in these directions die out ex- 
ponentially. Pure gauge fields on the lattice show a characteristic Lyapunov 
spectrum consisting of one third of each kind of exponents [5]. Assuming 
this general structure of the Lyapunov spectrum we investigate presently 
its magnitude only, namely the maximal value of the Lyapunov exponent, 

-^max- 

The general definition of the Lyapunov exponent is based on a distance 
measure d(t) in phase space, 

L := lim lim -In^T^- (5) 

t^oo d{o)^o t d{0) 

In case of conservative dynamics the sum of all Lyapunov exponents is zero 
according to Liouville's theorem, X^Lj = 0. We utilize the gauge invari- 
ant distance measure consisting of the local differences of energy densities 
between two 3d field configurations on the lattice: 

^•■=^T.P \t^Up-tvU'p\. (6) 

Here the symbol J2p stands for the sum over all Np plaquettes, so this 
distance is bound in the interval {0,2N) for the group SU(N). Up and 
Up are the familiar plaquette variables, constructed from the basic link 
variables Ux,i, 

Ux, = cxp (aAlT^) , (7) 



located on lattice links pointing from the position x = {xi,X2, xs) to x+aci. 
The generators of the group are = —igT'^/2 with t'^ being the Pauli 
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matrices in case of SU(2) and A"^ ^ is the vector potential. The elementary 
plaquette variable is constructed for a plaquette with a corner at x and 
lying in the ij-plane as Ux,ij = Ux,iUx+ijUl^j j. It is related to the 
magnetic field strength j^: 

Ux,ij = exp (^SijkaBl kT') . (8) 

The electric field strength E'^ ^ is related to the canonically conjugate mo- 
mentum Px i = Ux i via 



Ki = {t'Ux.uI) ■ (9) 



The Hamiltonian of the lattice gauge field system can be casted into the 
form 

hp,p) + i-hu,v) . 



(10) 



Here the scalar product stands for {A,B) = ^ii{AB'^). The staple variable 
y is a sum of triple products of elementary link variables closing a plaquette 
with the chosen link U . This way the Hamiltonian is formally written as 
a sum over link contributions and V plays the role of the classical force 
acting on the link variable U . 

Initial conditions chosen randomly with a given average magnetic en- 
ergy per plaquette have been investigated in past years [6]. In the present 
study we prepare the initial field configurations from a standard four di- 
mensional Euclidean Monte Carlo program on a 12^ x 4 lattice varying the 
inverse gauge coupling j3 [7] . We relate such four dimensional Euclidean lat- 
tice field configurations to Minkowskian momenta and fields for the three 
dimensional Hamiltonian simulation by identifying a fixed time slice of the 
four dimensional lattice. 



4. Chaos, confinement and continuum limit 

We start the presentation of our results with a characteristic example of the 
time evolution of the distance between initially adjacent configurations. An 
initial state prepared by a standard four dimensional Monte Carlo simula- 
tion is evolved according to the classical Hamiltonian dynamics in real time. 
Afterwards this initial state is rotated locally by group elements which are 
chosen randomly near to the unity. The time evolution of this slightly ro- 
tated configuration is then pursued and finally the distance between these 
two evolutions is calculated at the corresponding times. A typical exponen- 
tial rise of this distance followed by a saturation can be inspected in Fig. 2 
from an example of U(l) gauge theory in the confinement phase and in the 
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Figure 2. Exponentially diverging distance in real time of initially adjacent U(l) field 
configurations on a 12^ lattice prepared at /3 = 0.9 in the confinement phase (left) and 
at /? = 1.1 in the Coulomb phase (right). 




0.5 1 1.5 2 0.5 1 1.5 2 



Figure 3. Monopole density (left) and Lyapunov exponent (right) of the decomposed 
U(l) fields as a function of coupling. The Lyapunov exponent being the only observable 
with visible error bars. 



Coulomb phase. While the saturation is an artifact of the compact distance 
measure of the lattice, the exponential rise (the linear rise of the logarithm) 
can be used for the determination of the leading Lyapunov exponent. The 
left plot exhibits that in the confinement phase the original field and its 
monopole part have similar Lyapunov exponents whereas the photon part 
has a smaller L^^ax- The right plot in the Coulomb phase suggests that 
all slopes and consequently the Lyapunov exponents of all fields decrease 
substantially. 

We now turn to a comparison with the monopole density in U(l) quan- 
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Figure 4- Comparison of average maximal Lyapuriov exponents as a function of the 
scaled average energy per plaquette ag^E. The full U(l) theory shows an approximately 
quadratic behavior in the weak coupling regime which is more pronounced for the photon 
field. The monopole field indicates a linear relation. 

turn field configurations, 



Tlie left plot of Fig. 3 exhibits for a statistics of 100 independent configura- 
tions that Qmon decreases sharply between the strong and the weak coupling 
regime. It can be seen that the photon part from the decomposition carries 
also a few monopoles. 

The main result of the present study is the dependence of the leading 
Lyapunov exponent Lmax on the inverse coupling strength /?, displayed in 
the right plot of Fig. 3. As expected the strong coupling phase, where con- 
finement of static sources has been established many years ago by proving 
the area law behavior for large Wilson loops, is more chaotic. The transition 
reflects the critical coupling to the Coulomb phase. The right plot of Fig. 3 
shows that the monopole fields carry Lyapunov exponents of nearly the 
same size as the full U(l) fields. The photon fields yield a non-vanishing 
value in the confinement ascending toward /3 = for randomized fields 
which indicates that the decomposition (1) works perfectly for ideal con- 
figurations only. 

An interesting result concerning the continuum limit can be viewed from 
Fig. 4 which shows the energy dependence of the Lyapunov exponents for 
the U(l) theory and its components. One observes an approximately linear 
relation for the monopole part while a quadratic relation is suggested for 
the photon part in the weak coupling regime. From scaling arguments one 




X,fJ.U 



(11) 
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expects a functional relationship between the Lyapunov exponent and the 
energy [5, 8] 

L{a) DC a''-^E''{a), (12) 

with the exponent k being crucial for the continuum limit of the classical 
field theory. A value of A; < 1 leads to a divergent Lyapunov exponent, 
while k > 1 yields a vanishing L in the continuum. The case A; = 1 is 
special leading to a finite non-zero Lyapunov exponent. Our analysis of the 
scaling relation (12) gives evidence, that the classical compact U(l) lattice 
gauge theory and especially the photon field have A; ~ 2 and with L{a) — ^ 
a regular continuum theory. The monopole field signals fc ~ 1 and stays 
chaotic approaching the continuum. 

5. Conclusions 

We investigated the classical chaotic dynamics of U(l) lattice gauge field 
configurations prepared by quantum Monte Carlo simulation. The fields 
were decomposed into a photon and monopole part. The maximal Lya- 
punov exponent shows a pronounced transition as a function of the cou- 
pling strength indicating that on a finite lattice configurations in the strong 
coupling phase are substantially more chaotic than in the weak coupling 
regime. The computations give evidence that the Lyapunov exponents in 
the original U(l) field and in its monopole part are very similar. The sit- 
uation for the monopole density is analogous and serves as a consistency 
check of the decomposition. We conclude that classical chaos in field config- 
urations and the existence of monopoles are intrinsically connected to the 
confinement of a theory. 

We found evidence that the monopole fields stay chaotic in the contin- 
uum while the photon fields and the full U(l) fields possess a regular contin- 
uum theory. So far the monopoles are extracted from the Euclidean gauge 
fields. It will be interesting to compute their counterparts in Minkowski 
space and perform an analysis concerning monopole annihilation. 
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